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Abstract. We report ^'^^La, ^'^Fe, and ''^As Nuclear Magnetic Resonance (NMR) 
and Nuclear Quadrupole Resonance (NQR) measurements on powders of the new 
La0i_2;F2;FeAs superconductor for a;=0 and x=0.1 at temperatures up to 480 K, and 
compare our measured NQR spectra with LDA calculations. For all three nuclei in the 
a;=0.1 material, it is found that the local Knight shift increases monotonically with 
increasing temperature, and scales with the macroscopic susceptibility, suggesting a 
single magnetic degree of freedom. Surprisingly, the spin lattice relaxation rate for 
all nuclei also scale with one another, despite the fact that the form factors for each 
site sample different regions of q-space. This result suggests a lack of any q-space 
structure in the dynamical spin susceptibility that might be expected in the presence 
of antiferromagnetic correlations. Rather, our results are more compatible with simple 
quasi-particle scattering. Furthermore, we find that the increase in the electric field 
gradient at the As cannot be accounted for by LDA calculations, suggesting that 
structural changes, in particular the position of the As in the unit cell, dominate the 
NQR response. 
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1. Introduction 

The recent discovery [1] of superconductivity in the layered ferropnictides ROi-^F^^FeAs 
(R = rare earth) has raised great interest within the sohd state community. Not only 
does the transition temperature, Tc, reach a maximum at 55 K, but it is strongly 
dependent on the rare-earth ion and the pressure [21 [3], H]. Furthermore, the normal 
state properties exhibit some unusual features, which are reminiscent of the copper 
oxide high temperature superconductors (HTSC). In particular, there is a pseudogap- 
like decrease of the magnetic susceptibility at low temperatures [3 O [71 [8], and 3D 
antiferromagnetic order in the parent (undoped) compound LaOFeAs with T/v ~ 140 K. 
An important difference, however, is that the ferropnictides exhibit metallic properties, 
and clearly are not Mott insulators [9l [10] . Upon doping, the antiferromagnetic order 
is destroyed in both families, and a superconducting ground state emerges in the phase 
diagram. The proximity of superconductivity to an antiferromagnetic ground state and 
the appearance of the pseudogap features hint at the presence of magnetic correlations 
in these materials, as in high cuprates, that may play a critical role in the underlying 
physics of the superconductivity. 

In this article, we report Nuclear Magnetic Resonance (NMR), and Nuclear 
Quadrupole Resonance (NQR) investigations of the normal state of LaOo.gFg.iFeAs and 
LaOFeAs. Our results shed light on the role and importance of magnetic correlations 
in these compounds. In particular, we find no evidence for strong magnetic correlations 
in superconducting LaOo.gFo.iFeAs. NMR and NQR are well suited to probe magnetic 
correlations, since they are sensitive local probe techniques giving access to the intrinsic 
susceptibility, with a nucleus-dependent sensitivity to certain regions in q-space. In 
Sec. [2, we outline the experimental details of sample preparation and characterization, 
as well as NMR- and theory-related details. The NMR results on the electronic spin 
susceptibility are then presented and discussed in Sec. [3l Our Knight shift and spin 
lattice relaxation rate measurements of three different nuclei (^^As, ^^Fe, and ^^^La) 
scale with one another as a function of temperature. This result is surprising, as it 
suggests that all nuclei couple to the same spin degree of freedom, and that there is 
little or no g-space structure in the dynamical spin susceptibility. If spin fluctuations 
were present, with a correlation length greater than about one lattice spacing, then the 
^^As and ^^^La spin lattice relaxation rates would differ from that of the ^'^Fe, in contrast 
with our observations. We also measure the Knight shift of the ^^As up to 480 K, in 
order to investigate the behavior of the spin susceptibility at high temperature. We 
find that the shift increases monotonically up to 480 K, showing no sign of a peak as 
observed in pseudogap studies of the cuprates [H]. Finally, in Sec. [H we present density 
functional calculations of the spatial charge distribution and electric field gradient (EFG) 
at the ^^As site in LaOo.gFo.iFeAs and LaOFeAs and compare with our experimental 
observations. 
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2. Experimental details and theoretical methods 

2.1. Sample preparation and characterization 

Poly crystalline samples of LaOi_^F^FeAs with x=0.1 and a;=0 were prepared by 
standard methods and characterized by x-ray diffraction, resistivity and susceptibility 
measurements as described in [6l [12]. A value of ~ 26.0 K was extracted from 
these measurements for a fluorine doping of a;=0.1. A similar but ^^Fe-enriched sample 
showed a reduced Tc ~20 K from low-field SQUID measurement, with a reduced 
Meissner effect. The origin of the reduced is not yet clear. A change in the doping 
level could be excluded, as the ^^As NQR spectrum was unchanged. The undoped 
(non-superconducting) sample exhibits a structural phase transition from tetragonal to 
orthorhombic at Ts ~ 156 K followed by a spin density wave transition at T/v ~ 138 K 
[13]. For the NMR experiments, an oriented powder of the (unenriched) x = 0.1 sample 
was formed by grinding the material to a powder, mixing with Stycast 1266 epoxy, and 
curing in an external field of 9.2 T. 



2.2. NMR and NQR 

NMR is a powerful probe of the behavior of the electronic system in the ferropnictides 
because the nuclei interact with the electrons via quadrupolar and hyperfine interactions. 
The nuclear Hamiltonian is given by: 



H = -7,fil ■ (1 + K)Ho ■ 



3 6^ + 



1) 



2 

where 7„ and / are the gyromagnetic ratio and the spin of the nucleus, h is Planck's 
constant, Hq is the applied magnetic field, K is the NMR shift, and uq and rj the 
quadrupole frequency and the asymmetry of the electric field gradient tensor. The 
quadrupole moment, Q, of nuclei with spin J > 1/2 (La, As) interacts with the electric 
field gradient (EFG), which depends on the local charge symmetry. The quadrupolar 
interaction lifts the degeneracy of the / multiplet, giving rise to 2/ + 1 resonances in an 
applied field (NMR), and 7 — 1/2 resonances in zero field (NQR). 

The Knight shift K arises through the hyperfine interaction between the spins of 
the electrons and the spins of the nuclei. In the presence of an applied field (NMR), 
the static field of the polarized electrons creates an additional field at the nuclear sites, 
yielding a shift, K, of the resonance line with respect to the unshifted Larmor frequency. 
Furthermore, any time dependence of this hyperfine field (due to electronic spin and 
orbital moment fluctuations), gives rise to spin lattice relaxation. In general, K can be 
written as the sum of a temperature-dependent spin part. Kg, and a usually temperature- 
independent orbital part, Korb- Kg is proportional to the static susceptibility of the 
electrons at the Fermi level, XsiQ = 0): 

Kg = Ahf ■ Xs{q = 0) (2) 
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Table 1. NMR/NQR properties of the studied nuclei. The relative sensitivity for iron 
assumes 100% "Fe enrichment. 
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i39La 


7/2 


6.014 


0.2 


«100 


5.9-10-2 


"Fe 


1/2 


1.3757 





«2.2 


3.4-10-5 


75 As 


3/2 


7.2917 


0.314 


100 


2.5-10-2 



where Ahf is the hyperfine couphng, which depends on the nucleus through the 
availabihty of various couphng paths, i.e., on-site couphng or transferred couphng to 
neighbouring atoms. 

AU of the nuclei are NMR-active in LaOo.gFo.iFeAs; we have chosen to focus on 
^^^La, ^^As, and ^^Fe, whose properties are given in Tab. [H Since the La is located 
out of the FeAs plane, it is expected to be coupled only weakly to the electronic spin 
system. The EFG at the La site is particularly small, rendering NQR experiments 
difficult to carry out. On the other hand. As is located directly in the FeAs planes, 
and therefore should have a comparatively higher coupling to the electronic spins. 
Moreover, the combination of a large quadrupole moment and EFG allows for direct 
NQR experiments on the ^^As. The spin 1/2 nature of the ^^Fe is ideal for contrasting 
with the ^^As in order to distinguish quadrupolar from magnetic effects. However, the 
low natural abundance of the ^^Fe isotope makes it necessary to enrich the sample, 
using ^''Fe during the synthesis. Nevertheless, the sensitivity remains low compared to 
the other nuclei, making these experiments more difficult. 

When doing NMR on a powder sample, the crystallites are randomly oriented 
with respect to the applied magnetic field. This implies the random orientation of the 
magnetic and quadrupolar electric hyperfine tensors, giving rise to broad spectra and 
reduced signal intensities. Although EFG and magnetic shift tensors can be extracted 
from powder patterns in principle, we have chosen to exploit the anisotropic character of 
the magnetic susceptibility of these materials. By letting a mixture of powder and liquid 
epoxy cure while subjected to an external field, one obtains a powder with its crystallites 
oriented along the direction of highest susceptibility. In the case of LaOo.gFo.iFeAs the 
susceptibility is larger in the ab plane, i.e., along two directions. Therefore, we obtain ab 
oriented samples where the ab planes are parallel to the applied field, with a distribution 
of the resonance simpler (2D powder) than in the standard powder case. 

2.3. LDA calculations 

The band structure calculations were performed using the full-potential local-orbital 
minimum basis code FPLO (version 5.00-19) [11] within the local density approximation 
(LDA). In the scalar relativistic calculations the exchange and correlation potential 
of Perdew and Wang [T5] was employed. As basis set La (5s5p/6s6p5d-|-4f7s7p), Fe 
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(3s3p/4s4p3d+4d5s5p), As (4s4p3d+4d5s5p) and O (2s2p3d+3s3p) were chosen for 
semicore/valence+polarization states. The high lying states improve the basis which is 
especially important for the calculation of the electric field gradient (EFG) tensor with 
the components Vij = dV/dxidxj. The lower lying states were treated fully relativistic 
as core states. A well converged fc-mesh of 275 fc-points was used in the irreducible 
part of the Brillouin zone. LaFeAsO was calculated in space group 129 (P4/nmm) with 
the structural parameters as given in |9]. In order to investigate the influence of F 
substitution on the O site, the virtual crystal approximation (VGA) was applied and 
cross checked with the calculation of super cells [j. 

3. Spin susceptibility 

In this section, we present results pertaining to the electronic spin susceptibility of 
the FeAs layers, as probed by ^'^^La, ^^As, and ^''Fe nuclei. Through spectra and 
relaxation measurements, we probe the temperature dependence of the static and 
dynamic susceptibility. 

3.1. Uniform susceptibility 

3.1.1. Knight Shift Tensor at the Fe Accessing the temperature dependence of the 
susceptibility at q = was done through measurements of the NMR shift for the three 
studied nuclei, as derived from fitting of experimental spectra. A representative example 
is given on Fig. [H which shows ^^Fe NMR data for an unaligned doped LaOo.gFo.iFeAs 
powder. The data was obtained using the standard spin echo sequence with r=24 /is, 
while sweeping the field at a constant frequency of 12.7 MHz. These spectra correspond 
to the iron (1/2 -1/2) transition at several temperatures. For T=20 K (Tc < 20 
K at ifo = 9.2 T, see Sec. |2]), we can clearly resolve a high field tail in the spectrum. 
Since ^^Fe has / = 1/2, there are no quadrupolar contributions to this spectrum, and 
in fact this shape is typical of a powder distribution with an anisotropic Knight shift 
tensor. This interpretation is confirmed by a simulation of such a distribution (see Refs 
[inilTT]), shown as a solid fine in Fig. [H allowing us to extract the eigenvalues -ft'1/2/3 of 
the magnetic hyperfine tensor. We find Ki = K2 = 1.35(1) % and = 0.85(2) %, with 
the higher uncertainty on K3 due to the lack of definition of the end of the high-field 
tail {Ki/2 relate to the low-field peak). In light of the layered structure, it is reasonable 
to assume that the principal axes of the magnetic hyperfine tensor lie along the c axis 
and within the ab plane, with equal in-plane eigenvalues, and we assign K1/2 and K3 
as Kab and Kc- On increasing the temperature, we observe similar spectra, at least for 
the low-field peak (the poor sensitivity precludes comparison at higher field). The full 
lines show again powder simulations, with held at its value at T = 20 K so that the 

X In order to come close to the experimental F concentration of 10 %, a 4— fold super cell (doubled 
along a and b), with 8 formula units was calculated. Therefore the space group Pmm2 (#25) was 
chosen. Replacing one O by F yields a composition of LaFeAsO0.875F0.125- 
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high-field decrease remains the same. The low-field peak, which is solely determined by 
Kab, shows a weak temperature dependence (see Fig. [2]). If we assign the temperature- 
independent component to the orbital shift (as is customarily done), we find a quite 
small (less than 0.1%) spin component of the NMR shift along ab, with most of the 
shift being of orbital origin {K^^^ ^ 1.4%). Note however the much smaller value of Kc 
compared to Kab- In a similar ^^Fe study on LaFeAsOo.7, Terasaki et al [18] measure 
an even bigger anisotropy at T=30 K, with a shift of roughly 0.5% along c compared to 
again 1.4% in the plane. One possibility would be that this strong difference between 
the c and ab directions originates from a large difference in spin shift, well exceeding 
an order of magnitude, with an isotropic orbital shift. This would imply a very large 
anisotropy of the iron hyperfine coupling, since SQUID measurements rule out any such 
large anisotropy of the susceptibility. Having an isotropic orbital shift in spite of the 
doping would also likely imply at minimum significant mixing of the 3d orbitals. While 
we cannot rule out that scenario, a more straightforward explanation would be that 
there is a large orbital shift anisotropy (implying K^^.^ f» 0.8%), refiecting the fact that 
the iron valence shell is not closed. This may yield indications on the average valence 
of the iron, thus on the doping, as was suggested by Mukhamedshin et al [I9] in the 
study of sodium cobaltates. 

3.1.2. Scaling and Evidence for a Single Component The extracted temperature 
dependence of the ^^Fe shift Kab is presented on Fig. [21 Also shown are the ^^As 
(extended to high temperatures from previous work [6]) and ^^^La Knight shift in the ab 
directions, Kab, as measured in oriented LaOo.gFo.iFeAs powder (2D powder) with a field 
ifo=7.0494 T applied along the a, b directions. The resonance frequencies of the central 
transitions of the As and La are given to second order by / = Hq{1+ K ab) + 'iu'Q / IQ-^ Hq, 
with / the frequency of the peak of the 2D powder spectrum. For ^^As, we have 
independently measured vq = 11.00(5) MHz (see Fig. [Hand Ref. [6]), and therefore 
we can extract Kab from the spectrum. For the La, we measured the position of the 
satellites in a full spectrum, and found that uq = 1.15(5) MHz. The macroscopic powder 
susceptibility measured in an applied field of 5 T is shown as a solid line. By proper 
scaling, the four data sets can be made to overlap in the paramagnetic region. The 
legitimacy of such a procedure is based on the fact that each shift can be written as the 
sum of a temperature- independent orbital shift Korb, plus a spin shift Kg. Likewise, the 
static susceptibility can be written as the sum of a diamagnetic term Xdia plus a Van 
Vleck-like term Xvv, both expected to be temperature- independent, and a temperature- 
dependent spin term Xs,macro- Therefore the ratio of the scales of the NMR shifts to the 
scale of Xpowder rcficcts the strength of the hyperfine couphng of each nucleus to the spin 
susceptibility. It is remarkable that such general scaling can be obtained, and suggests 
that the three nuclei probe the same component of spin susceptibility. 

A priori this result is surprising, as there are several bands crossing the Fermi 
surface in the ferropnictides, and one might expect each band to supply a different 
contribution to the spin susceptibility with different hyperfine couplings to different 
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Figure 1. Field sweeped NMR spectra of tiie ^ Fe for different temperatures measured 
at an applied frequency of 12.7 MHz. The full lines are simulations as described in the 
text. 
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Figure 2. Knight shift of ^^As (black squares), ^^Fe (green down-pointing triangles) 
and ^^^La (blue triangles), and the macroscopic susceptibility Xpowder (solid red line) 
versus temperature with different vertical scales and origins. Note the reversed scale 
for 57Fe. 



bands. Such is the case for the oxygen in Sr2Ru04, which simultaneously couples to 
multiple bands with different temperature dependent susceptibilities [20]. The Mila- 
Rice-Shastry picture in the cuprates, however, captures much of the relevant physics in 
terms of a single spin component [21], [221 ESI El]- Our results suggest that if there are 
different hyperfine couplings to multiple orbitals in the ferropnictides, the spin response 
of each of it is nearly identical. This observations may reflect the itinerant, rather than 
localized character of the Fe 3d electrons. 

The shift data in Fig. 2 clearly show a strong decrease of the spin susceptibihty 
with decreasing temperature, in agreement with several previous reports [SlEllH]. This 
suppression of low-energy spin excitations is similar to the behavior of the cuprates, and 
hence has been ascribed to the existence of a pseudogap in this system. In the cuprates, 
the spin susceptibility reaches a maximum at a temperature T* that is doping dependent. 
Since our Knight shift measurements indicate an increasing susceptibility with increasing 
temperature, we sought to find if this trend continues to higher temperature. Our 
^^As measurements up to 480 K show no signature of a pseudogap peak, although 
the data hint at a leveling out by 500 K. It is unclear whether points at even higher 
temperatures could be gathered, as was for instance fruitful for YBa2Cu408 [H], since 
SQUID measurements suggest a degradation of the sample. Note that the scaling of the 
NMR shift with Xpowder remains good down to T^, which indicates the high quality of the 
samples, and the absence of any signature of a Curie contribution due to paramagnetic 
impurity spins, either intrinsic or belonging to a spurious phase. 
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Table 2. Hyperfine couplings and orbital shifts for ^^^La, ^'^Fe, and "^^As in the 
LaOo.gFo.iFeAs compound. The results for ^^As are from previous work [5]. 





i39La 


57pg 


75 As 


Kt, (%) 


4.3(8) 
0.12(1) 


-5.7(14) 
1.36(1) 


25(3) 
-0.03(4) 



3.1.3. Hyperfine couplings Table [2] summarizes the hyperfine couphngs and 
temperature-independent shifts extracted by plotting the measured shifts versus the 
bulk susceptibility, using Eq. [2l Since the anisotropic components of x ^^re unavailable 
in our aligned powder, we compare with the powder susceptibility in this determination, 
ignoring the modest anisotropy in temperature dependence between the c and ah 
directions [25]. We do not have an estimate for the non-spin component Xvv + Xdia 
but it is between an overly-cautious lower bound of and an upper bound of roughly 
2 ■ 10~^ emu/mol, otherwise the spin contribution would have to become negative below 
a certain temperature. This gives for the temperature-independent fraction of the shift 
^^^K^\ = 0.12(1)% and K'^\ = 1.36(1)%. 

The largest hyperfine coupling is to the As, which may be due in part to the fact 
that there are four nearest neighbor iron atoms to each As. The fact that ^^'^Ahf is 
roughly six times lower than '^^A^f is not surprising given that lanthanum is outside 
of the FeAs layers. In the case of iron, the negative hyperfine coupling suggests that 
the dominant hyperfine coupling is via a core polarization mechanism, but the small 
magnitude of the Fe hyperfine coupling is surprising. A priori., one would expect that 
iron would have the strongest coupling to the electronic properties, in light of theoretical 
predictions [26l [271 [2H] indicating the highly predominant iron "id character of the bands 
at the Fermi level. One explanation is that ^^Ahf = A^p + where Acp (the core 
polarization contribution) is large and negative, whereas is large and positive, so 
the net coupling remains small. 

From the measured lanthanum coupling, it is possible to estimate the value of 
the iron moment m^g in the magnetically- frozen SDW state of the undoped LaOFeAs 
material. As there are no major structural changes between the 10% fluorine- doped 
material and the undoped material in its low-temperature state (aside from the slight 
orthorhombic distortion as a prelude to the magnetic transition) it is reasonable to 
assume that ^^^A^f remains doping independent. In this case, using the internal field of 
Hint{La) = 2.5 kOe measured by Nakai et al at the lanthanum site in the SDW state 
[7j and the relation Hint{La) = ^^^Ahfmpe, we estimate mpe = 0.58(9)/xb. While this 
is larger than the values obtained through Mossbauer [10] and neutron scattering [9] 
measurements, which gave respectively mpe =0.25 fiB and mpe =0.36 fiB, this result 
still points clearly at a largely itinerant situation in the SDW state. 
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Figure 3. Knight shift of ^'^^La (a) and of ^^Fc (b) versus the macroscopic 
susceptibihty Xpowder- The black hnes are linear fits. The orbital shifts (see text) 
are shown next to the left vertical axes. 



3.2. Wavevector dependence of the susceptibility 

As the superconductivity and frozen magnetism regions are close to each other in the 
phase diagram, it is natural to consider that spin fluctuations may play a role in the 
formation of Cooper pairs. NMR measurements of the spin lattice relaxation T^^ sample 
the low energy spin fluctuations via the dynamic susceptibility x" 

r-oc%^limV|^(,)r^l(«l^, (3) 

where /c^ is Boltzmann's constant, 7e/n the gyromagnetic ratios of the electron and the 
probed nucleus, Uq the Larmor frequency, and A{q) the hyperfine form factor of the 
probed nucleus. Particular attention should be paid to the g-dependence of the latter, 
as filtering effects may occur such as in the superconducting cuprates [30], wherein the 
oxygens in the Cu02 layers are insensitive to antiferromagnetic fluctuations, whereas 
the copper nuclei are. Here, the fact that ^''Fe appears to be a particularly poor probe 
of the uniform {q=0) susceptibility compared to ''^As and ^^^La, taking into account 
the crystallographic positions, suggests that hyperfine filtering effects might be at play. 
While the numerous superexchange paths make the analysis difficult for fluorine and 
lanthanum, Terasaki et al [18] show indeed that '^^Ahf{q) is well-developed at q = and 
tends to vanish towards q = (7r/a,7r/a), while ^"^Ahfiq) exhibits an opposite behaviour, 
tending to low values around q = 0. 

The temperature dependence of (TiT)~^ for each of the Fe, As and La nuclei is 
shown on Fig. HI including ^^F data from Ahilan et al [5] for comparison. For ^^^La, 
^^Fe, and ''^As, Ti along the ah directions is used, while for fluorine it is T^*'^°. The data 
are plotted with different axes, in order to highlight the similar temperature dependence 
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Figure 4. The temperature dependence of the ^"^Fe, ''^As, ^^^La, and ^''F (rir)"^ 
The ^^F data is reproduced from Ahilan et al [5]. For -'^^^La, ^^Fe, and ^^As, Ti along 
the ab directions is used, while for fluorine it is Tl'^° . 



for all four nuclei. The strength of the relaxation correlates with the distance to the 
iron plane, with the lowest values for the nuclei outside the FeAs planes. The spin 
lattice relaxation rate is largest for the Fe, whereas the spin shift at the Fe is relatively 
small compared to the other sites. In other words, ''^(TiT7^)/^''(TiT7^) ?«20-30 while 
{^'^AhfY{q = 0)/(''^A/j/)^(q = 0) «i0.05. There are two possible explanations for this 
discrepancy. Either there is a strong g-space dependence to x"{q,uj), or there are 
multiple hyperfine coupling channels {Acp, ^4^) between the Fe nuclear moments and 
the same degree (s) of spin freedom. Given the fact that the spin lattice relaxation rate 
for all the nuclei exhibit roughly the same temperature dependence as seen for the NMR 
shifts, any strong q dependence to the dynamic susceptibility seems unlikely. In other 
words, while the susceptibility is temperature-dependent, the decrease occurs similarly 
across q-space. While some fluctuations in certain q regions cannot be ruled out in the 
absence of a refined analysis of hyperfine filtering effects, a simple explanation of our 
data would be that the relaxation comes mostly from quasi-particle scattering. In this 
case, the core-polarization and the diamagnetic contributions to the Fe relaxation rate 
add as the sum of the squares rather than a direct sum, as discussed in [31] . We note 
that there are no signatures of a pseudogap peak in the (TiT)~^ data, as observed in 
the cuprates [211 122] in either the Fe or As data. 
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4. Spatial charge distribution 

As a complementary study, we present in this section a theoretical analysis of the issue 
of the spatial charge distribution, based on our ''^As NQR measurements. The EFG 
observed at the As site in the ferropnictides varies dramatically from one compound 
to the next, in stark contrast to the cuprates, where the EFGs are slightly doping 
dependent, but exhibit little variation among different families. To address this, we 
have measured the doping dependence of the NQR spectrum, and compared with LDA 
calculations. 

4.1. NQR results 

In the ferropnictides, the ^^As nucleus quadrupolar moments interact with the electric 
field gradient (EFG) tensor giving rise to a resonance at frequency: 



where r/ is the asymmetry of the EFG tensor. The EFG has the symmetry of the local 
atomic position, and depends on the local electronic density. We present on Fig. [5] 
the ^^As NQR spectrum at room temperature for the doped LaOo.gFo.iFeAs material 
studied here above, as well as for the undoped parent compound LaOFeAs. For both 
dopings, a well defined line is observed, meaning that in each sample the EFG is the 
same at all As sites. This is in agreement with the single As crystallographic site 
and indicates spatially homogeneous doping. In the doped sample, the line is however 
significantly broader (full width at half maximum of 0.97(6) MHz) than in the undoped 
case (FWHM=0.22(1) kHz). This could be explained by limited inhomogeneities of the 
fluorine concentration in the material, distributing somewhat the EFG, or even simply 
by the fact that As ions are at varying distances from the fluorine. We find that uq 
increases from 9.48(1) MHz to 11.00(5) MHz on doping, which translates accordingly in 
a 16% increase of Vzz- This trend is in good agreement with measurements by Mukuda 
et al [33] on oxygen-deficient compounds. 

4-2. Comparison with theory 

For the undoped compound the measured EFG is obtained by inserting the ^^As 
quadrupole moment [3l] Q = (3.14 ± 0.06) b in Eq. (jlj), using ?7=0, which yields 
IKT^I = (2.50 ± 0.05) ■ 10^^ V/m^. Using the 175 K lattice parameters and atomic 
positions as given in Ref. [9] we obtain a fair agreement for the calculated EFG: 
ycaic _ _3 . iq2i v/m^. Like for the Fe magnetic moment, which in the calculations 
shows a strong dependence on the As z position [35l [36] , we observe also for the EFG a 
strong As z dependence (inset to Fig. [6]). The calculated EFG {V^^f'' = -2.67- lO^^ V/m^) 
is much closer to the measured value when As is shifted along the negative z-direction 
to z = 0.6438, where the energy has a minimum (in inset to Fig. [6] marked by an arrow) 
and the structure has a shorter Fe-As distance of 2.3748 A. 




(4) 
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Figure 5. Room temperature ^^As NQR spectra in tlie doped LaOo.gFg.iFeAs 
material and in the undoped parent compound LaOFcAs. Solid lines indicate gaussian 
fits, with parameters given in the text. 
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As z displacement [A] 

qL , 1 , 1 , 1 , 1 

0.05 0.1 0.15 0.2 

doped electrons 

Figure 6. Calculated Vzz obtained from the virtual crystal approximation using the 
experimental As z = 0.6507 position (empty squares), the optimized As z = 0.6438 
position (filled diamonds) and the LaFeAsO0.92F0.08 structure at 175 K (black circle). 
The red filled circles show the EFGs from the super cell calculation for the optimized 
As z position. The measured EFGs for the pure and the 10 % F-doped compound 
are shown by error bars. Inset: Dependence of ^^As V^z on the As z position. The 
energetically optimized As z position is marked by an arrow. The experimental Vzz is 
represented by the shaded bar. 

The EFGs of the doped compounds were calculated with the virtual crystal 
approximation. The validity of the VGA was confirmed by super cell calculations. 
Due to the super cell construction, there are two different Wyckoff positions for As and 
hence two different EFGs, whereof one is lying on top of the VGA curve and the other 
one very close to it, see full red circles in Fig. [61 First, we consider solely the effect 
of electron doping. Therefore, we keep the structural parameters fixed for different 
levels of doping. In Fig. [6] two such VGA curves are shown. When the experimentally 
determined As z = 0.6507 position is used, the calculated and measured EFG values 
for 10 % doping agree quite well. Also S. Lebegue et al found good agreement for the 
10 % doped compound using the WIEN2k code [37]. The VGA curve with the optimized 
As z = 0.6515 position is shifted in the direction of smaller \Vzz\- Now, we investigate 
the structural change on top of the doping by calculating the EFG within VGA for the 
175 K data of LaFeAsO0.92F0.08 as given in Ref. [9]. This has only a minor effect on 
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the EFG, as it can be seen by the black circle in Fig. 0, which lies very close to the 
experimental VGA curve. We stress at this point, that the effect of electron doping on 
the EFG is much smaller than the influence of the As z position as can be clearly seen 
by comparing Fig. [6] with the inset of Fig. [61 

Our calculations result in a decrease in \Vzz\ upon electron doping for 
LaFeAsOo.gFo.i, although in our experiments an increase is observed. This is not pointed 
out by S. Lebegue et al [37], although they obtain the same discrepancy for the trend 
in the Vzz calculation. Further studies are required to investigate the connection with 
intrinsic changes in the electronic structure. 

For the NdFeAsO system, note that there is a better agreement between the 
experimentally determined and the calculated EFG for the doped and for the undoped 
structure [38]. In that study, some of us (K. K., H. R) also showed, that the 4/ electrons 
have only a minor influence on the EFG, whereas the structure (chemical pressure) 
influences the EFG to a higher degree. 

5. Conclusion 

Using simultaneous ^^^La, ^^Fe, and ^^As NMR and NQR measurements, we have 
investigated the electronic properties of the LaOi^^jFajFeAs compound with x=0 and 
x=0.1. ''^As NQR measurements show a sizable evolution of the electric field gradient 
with doping that cannot be explained by LDA calculations, although the measured and 
calculated EFG are in reasonable agreement for the undoped parent compound. While 
a high sensitivity to the As z position is observed, the electronic origin of the difference 
to the experimental spectra remains unclear. The systematic scaling of the NMR shifts 
with the macroscopic susceptibility in the paramagnetic state suggests the presence of a 
single spin degree of freedom. Extending ^^As shift measurements up to 480 K reveals a 
continuous increase of the spin susceptibility with no sign of a peak, suggesting that the 
observed pseudogap behavior persists up to at least this high. We find that the hyperfine 
coupling to the iron is unexpectedly small, but speculate that there are in fact two 
components to the Fe hyperfine coupling with different signs. Gonsequently the Knight 
shift at the Fe is rather small, but the spin lattice relaxation rate is large. Furthermore, 
we find that the spin lattice relaxation rates for all four nuclei appear to scale with one 
another, providing further support for a single spin degree of freedom, and suggesting 
the absence of any significant structure in the dynamic susceptibility in qr-space. This 
result is surprising: in contrast with the cuprates, there appear to be very little spin 
fluctuations at low energies present in these superconducting samples. Detailed studies 
of the NMR relaxation and its doping dependence, as well as quantitative comparisons to 
inelastic neutron scattering data, should help bring further insight on the issue of the role 
of spin fluctuations in these materials. The origin of the strong temperature dependence 
of the susceptibility, as well as the presence of unconventional superconductivity in the 
absence of significant spin fluctuations remain open questions, and indicate that physics 
of the iron pnictides are very different from that of the cuprates. 
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